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Abstract 

We show that the coupled mode equations for the stationary propagation of upper-hybrid and 
magnetoacoustic waves in magnetized electron-ion plasmas with negative group dispersion can 
be exactly derived from the generalized Henon-Heiles Hamiltonian. The parameter regimes for 
the integrable cases of the coupled mode equations have been explicitly obtained. For positive 
group dispersion of the upper-hybrid waves, the relevant governing equations lead to a novel 
Hamiltonian where the kinetic energy is not positive definite. 
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1. Introduction 



Amplitude modulated upper hybrid waves in a magnetized electron-ion plasma are known to 
be governed by a Schrodinger-like equation wherein the potential is given in terms of the associated 
low-frequency number density perturbations [1-3]. For small amplitudes, the latter are governed by 
a linear wave equation driven by the ponderomotive force of the high-frequency upper-hybrid waves. 
On the other hand, for finite amplitudes, the density perturbations are governed by a driven nonlinear 
Boussinesq equation [4] which is coupled to the Schrodinger equation. For uni-directional propagation, 
the (driven) Boussinesq equation reduces to the well-known (driven) Korteweg-de Vries (K-dV) equa- 
tion. For stationary propagation of the coupled waves, the time-dependent Schrodinger-Boussinesq 
(or, K-dV) equations give rise to a coupled system of nonlinear ordinary differential equations which 
can be derived from a Hamiltonian. While special analytical solutions valid in some specific regions of 
the allowed parameter space have been obtained [5], the question of the complete integrability of the 
associated Hamiltonian for arbitrary boundary conditions has not yet been answered. 

We show, in this Letter, that the coupled mode equations for the upper-hybrid and magnetoa- 
coustic waves with negative group dispersion can be exactly reduced to the generalized Henon-Heiles 
equations which are extensively studied in the field of Nonlinear Dynamics [6-8] . We thereby obtain 
explicitly the parameter regimes for the integrability of the associated Hamiltonian. For the case of 
positive group dispersion, the coupled mode equations are derivable from a novel kind of Hamiltonian 
having indefinite kinetic energy. 



2. Governing Equations 

We consider the one-dimensional propagation along x-direction of the high-frequency upper- 
hybrid waves in a homogeneous, magnetized electron-ion plasma with the external magnetic field (B Q ) 
along the ^-direction, that is, B = B z. For normal modes, the upper-hybrid wave frequency (ujo) 
and the wavenumber (k ) are related by the linear dispersion relation [9], 

%,) 2 k 2 v 2 
^0 ~ 4 "e0 

where u pe Q = (47rn e 2 /m e ) 1 / 2 is the electron plasma frequency, S7 e0 = eB /m e c is the electron gyro- 
frequency, v te = (Tg/mg) 1 / 2 is the electron thermal speed, to H0 = (w 2 e0 + ^eo) 1 ^ 2 ^ s ^ ne upper-hybrid 
frequency, and all the other symbols have their usual meanings [4]. In the long wavelength limit, Eq. 
(1) can be approximated by 
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where 



uj = uj h0 + i D Q k 2 , (2) 



_ d 2 U 3UpeO V te ,os 

u ® = "hit - ~. — t~72 ' \ 6 ) 



dk 2 uj H0 {u 2 e0 - 3tt 2 



c.OJ 



denotes the group dispersion coefficient for the upper-hybrid waves. Clearly, the latter have positive 
(negative) dispersion for plasma parameters such that u 2 e0 > 3f2g (oo 2 eQ < 3VL 2 ). 

For nonlinear propagations, the slowly varying complex amplitude E(x,t) of the upper-hybrid 
wave electric field is governed by a Schrodinger equation of the form [2,4], 

. (dE Tr dE\ D d 2 E 

where, iV = 5n e /rio is the normalized low-frequency density perturbation, V g = duo/dk = k D 
denotes the group velocity and fx = ^(uj 2 e0 + 2VL 2 )/(uo 2 e0 + Qg ). 

For finite amplitudes, the low-frequency density perturbation (JV) is governed by a Boussinesq or 
a K-dV type of nonlinear equation which is driven by the ponderomotive force of the high-frequency 
upper-hybrid waves. The driven Boussinesq or the K-dV equation can be derived from the low- 
frequency fluid equations (for the electrons and the ions) which are coupled to the Maxwell equations. 
Omitting the details of the derivation [4] , we write below the driven Boussinesq equation in the form, 



2 JL-V 2 — -R 2 — - ^(N 2 \= 2 ^ ( |E|2 ^ 



where, V M = (V 2 + C 2 ) 1 ^ 2 is the magnetoacoustic speed, V A = (B 2 /ATrn mi ) x / 2 is the Alfven speed, 
C s = (T e /rrii) l l 2 is the ion-acoustic speed, 9 = cV A /ujp e o, a 2 = (3V 2 + 2C 2 )/2 ) and rj = uj HO C s /uj pe o. 
Note that in the linear limit for normal modes, Eq. (5) gives, 

u 2 = V 2 k 2 - 9 2 k' , (6) 



which is the linear dispersion relation for the magnetoacoustic modes in the long wavelength regime. 
Equation (5) describes the bi-directional propagation of the nonlinear magnetoacoustic waves driven 
by the upper-hybrid ponderomotive force. 

On the other hand, for uni-directional propagation, Eq. (5) can be reduced to the form, 
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~dt + M ~dx~ + 2V^ ~dx T + ~dx " ~2V^dx~ \l&irn T e ) ' ^ 

which is the driven K-dV equation. For normal modes, Eq. (7) yields the dispersion relation, u = 
V M k — 9 2 k 3 /2V M which follows also from Eq. (6) in the small wavenumber limit. 

For stationary propagation, the wave fields are expressed in form, 

E(x,t) = E(Z)exp[i{X(x)+T(t)}] , (8.1) 

N(x,t)=N(0, (8.2) 

where £ = x — Mt represents the coordinate in the stationary frame whose speed is determined by the 
free parameter M; the functions X(x) and T(t) are introduced in Eq. (8.1) to account for the possible 
shifts in the wavenumber as well as in the wave frequency due to the nonlinear interactions. Using 
Eqs. (8.1) and (8.2) in Eq. (4), we obtain, 

D ^ = \E + b 2 NE , (9) 

where A = 25 + (M 2 — Vg)/D is the nonlinear shift parameter, 5 = dT/dt denotes the shift in the 
wave frequency and b 2 = 2 huj h q. 

The stationary governing equation for the density perturbations (JV) is obtained from Eq. (5) or 

(7) as, 

where, 

!M 2 — V 2 , for the driven Boussinesq Eq. (5), 

(11) 
2V M {M - V M ) , for the driven K-dV Eq. (7). 

Equations (9) and (10) are the relevant coupled mode equations for the stationary propagation 
of upper-hybrid and magnetoacoustic waves. We have presented elsewhere [5,10] different classes of 
exact analytical solutions which either use special boundary conditions or are valid in limited regions 
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of the allowed parameter space. In the next section, we reduce these equations to the generalized 
Henon-Heiles form and thereby determine the integrable parameter regimes. 

3. Reduction to Henon-Heiles equations 

In order to reduce the coupled mode equations (9) and (10) to the generalized Henon-Heiles form, 
we note that the parameters A, D and / can have either signs. We normalize N with respect to 
-2D D/b 2 and E 2 with respect to lQirn T e (2\D \D/b 2 ), and obtain, 

dPE - X F 2T]2 FN (12) 
d 2 N f Ar 2a 2 rf \D Q \ Ar2 r) 2 ~ 

where p — — 1 for negative dispersion (D < 0) and p — +1 for positive dispersion (D > 0) of the 
upper-hybrid waves. We consider below these two cases separately. 

(A) Negative dispersion (D < 0) 

For upper-hybrid waves with negative group dispersion (u 2 e0 < 3f2g ), the coupled mode equa- 
tions (12) and (13) become, 



d 2 N 



d 2 E 

■ —AE — 2DEN , (14) 



-BN - CN 2 - DE 2 . (15) 



de 

where, 

a= \d;v ^' c = —^> D = ¥- (16) 

Equations (14) and (15) can be derived from the Hamiltonian, 

H + = l - (n 2 E + U 2 N ) + l - (AE 2 + BN 2 ) + QcW 3 + DNE 2 ^ , (17) 

where U E = dE/d^ and = dN/dS, are, respectively, the "canonical momenta" conjugate to E and 
N. 
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Clearly, by treating E and N as the spatial coordinates and £ as the temporal coordinate, Eq. 
(17) may be considered as the Hamiltonian for the two-dimensional motion of a pseudo-particle of unit 
mass. In fact, H + is identically the same as the generalized Henon-Heiles Hamiltonian [7,11] which has 
been extensively studied in the field of Nonlinear Dynamics. Since the stationary coordinate (£) does 
not explicitly appear in Eq. (17), the associated potential is conservative and hence the Hamiltonian 
H + is an integral of motion. For the two-dimensional motion, the system is completely integrable 
provided there exists the second integral of motion which is in involution with the Hamiltonian [12]. 
On the other hand, in recent years, the so-called "Painleve Analysis" has been extensively used to 
obtain the parameter regimes wherein low-dimensional Hamiltonian systems may be completely inte- 
grable [7,8]. In particular, it is well-known that the Hamiltonian (17) is completely integrable for the 



following three sets of parameter values : 

(a) A = B, C = D (18.1) 

(b) 16A = B, C = 16 D (18.2) 

(c) arbitrary A and B, C = 6D. (18.3) 



The associated second integrals of motion for the above parameters have been summarized in Ref. 
[8]. Furthermore, there are indications from general considerations that these are possibly the only 
integrable cases of the generalized Henon-Heiles Hamiltonian [11]. Thus, a necessary condition for the 
integrability of Eqs. (14) and (15) seems to be that the nonlinear terms should have the same signs. 

To obtain explicitly the various plasma parameters for integrability, we shall first consider the 
case when Eq. (10) together with (11) corresponds to the driven Boussinesq equation (5). In terms 
of the dimensionless parameters defined by a = u pe o/fl e o, P = (C S /V A ) 2 and 7 = v te /c, the above 
relations can be written in the form, 

3aV(3 + 2£) = z/(2 + a 2 )(3 - a 2 ) , (19) 

and 

3«V(1 + /?)(! - M 2 ) = uA(l + a 2 )(3 - a 2 ) , (20) 

where A = X/uj h0 is the normalized nonlinear shift parameter, M is the Mach number normalized with 
respect to V M and v takes values 1, 6, or 16. Note that the parameter f3 is essentially the usual plasma 
beta, that is, the ratio of the thermal pressure to the magnetic field pressure. For integrability of the 
coupled Eqs. (14) and (15), both the conditions (19) and (20) should be simultaneously satisfied for 
the case when v — 1 or 16, whereas only the condition (19) needs to be satisfied when v — 6. 
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Defining T 2 



7 2 /z/, Eq. (19) can be written in the form, 



(2 + q 2 )(3-a 2 ) 
3a 4 (3 + 2(3) 



(21) 



Note that T 2 remains positive definite since a 2 < 3 is satisfied for negative group dispersion. Substi- 
tuting for 7 2 from Eq. (19) into Eq. (20), we get, 



Figure (1) shows a plot of T 2 as a function of a 2 from Eq. (21) for different values of (3. The value 
(3 = corresponds to the cold plasma case whereas large values of (3 correspond to weakly magnetized 
plasmas. For v — 6, the system of equations (14) and (15) is completely integrable for parameters 
given by Figure 1 and which satisfy T 2 = (3/a 2 v since (3 = a 2 7 2 . Note that for v — 6, any arbitrary 
values of M and A are admissible. On the other hand, Figure (2) gives a plot of A as a function of a 2 
from Eq. (22) for different values of (3. For v — 1 and 16, the corresponding values of the parameter 
7 for integrability are given by Figure 1 whereas the parameters M and A are no longer arbitrary but 
are related by Eq. (22). Since the right-hand side of Eq. (22) is positive definite always, it follows 
that for positive (negative) values of the frequency shift parameter A, the governing equations (14) 
and (15) are integrable for sub-magnetoacoustic, that is, for M < 1 (super-magnetoacoustic, M > 1) 
values of the Mach number M. 

The parameter values for integrability for the driven K-dV case in Eqs. (10) and (11) can 
similarly be obtained. In fact, it follows by inspection that Eqs. (21) and (22) hold good in this case 
also provided the factor (1 — M 2 ) in the latter is replaced by 2(1 — M). The Mach number (M) regimes 
for integrability are, therefore, qualitatively the same as in the driven Boussinesq case. 

(B) Positive dispersion (D > 0) 

For u 2 e0 > 3fi 2 , that is, for a 2 > 3, the upper-hybrid waves have positive group dispersion and 
the coupled mode equations (12) and (13) become, 



A = 



1 — M 2 
A 



(l + a 2 )(3 + 2(3) 
(2 + a 2 )(l+p) ' 



(22) 



d 2 E 



— AE — 2D EN , 



(23) 



d 2 N 



-BN + CN 2 + DE 2 , 



(24) 



de 



where the coefficients A, B, C and D are defined, as earlier, by Eqs. (16). The coupled equations (23) 
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and (24) are derivable from the Hamiltonian, 

H- = \ (n 2 B -K)~\ (AE 2 + BN 2 ) + QcW 3 + DNE 2 ^ , (25) 

where the canonical momenta, for the present case, are given by U E = dE/d£ and IT^ = —dN/dC,. As 
earlier, the Hamiltonian if_ is an integral of motion corresponding to Eqs. (23) and (24). 

The Hamiltonians H+ and if_ given, respectively, by Eqs. (17) and (25) differ from each other by 
sign changes for the terms containing U 2 N , E 2 and iV 2 . The sign changes for the terms containing E 2 
and iV 2 are trivial and can, in fact, be absorbed by redefining the coefficients A and B. However, the 
sign change for the quadratic term in is indeed significant. For, unlike the case of the Hamiltonian 
H + , the "kinetic energy" term in the Hamiltonian H_ is not positive definite. The latter is in contrast 
to the usual Hamiltonian systems of classical dynamics where the kinetic energy is positive definite 
always. An important consequence of this difference is that systems with indefinite kinetic energy 
need not necessarily have bounded motions around those points where the potential energy has a 
minimum. In fact, for such systems, both the canonical momenta can simultaneously increase or 
decrease but still keeping the Hamiltonian an integral of motion. Furthermore, the usual stability 
theorems developed in classical dynamics may not be directly applicable to such cases. It is therefore 
expected that the qualitative nature of the solutions in both the cases should be different. Hamiltonians 
with indefinite kinetic energy [13] and having different kinds of potential functions are known to arise 
in many problems dealing with the nonlinear evolution of the modulational instability of an high- 
frequency wave coupled to a suitable low-frequency wave. We have reported earlier [5,10] some special 
classes of exact analytical solutions of the coupled equations (23) and (24). However, the question 
of the complete integrability of such Hamiltonians, in general, and that of H_ given by Eq. (25), in 
particular, by either the standard Painleve analysis or otherwise is still open. 
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4. Conclusions 



To conclude, we have shown that the stationary governing equations for modulated upper-hybrid 
waves coupled to magnetoacoustic waves in a magnetized plasma with negative group dispersion are 
derivable from the generalized Henon-Heiles Hamiltonian of nonlinear dynamics. The parameter 
regimes for the integrability of the associated Hamiltonian have been explicitly obtained. For the 
case of upper-hybrid waves with negative group dispersion, the equations give rise to a novel type of 
Hamiltonian with indefinite kinetic energy. The Painleve analysis as well as the integrability of the 
latter remains to be investigated. 

The results of the present investigation should have a bearing on the possible parameter regimes for 
the occurrence of upper-hybrid turbulence in magnetized plasmas [3,14,15]. For example, magnetized 
plasma turbulence has been invoked as a possible source for the narrow-band, non-thermal continuum 
radiation observed in the Earth's magnetosphere [16]. In the present work, we have explicitly obtained 
those parameter regimes where the stationary wave fields are non-chaotic and hence can give rise to 
soliton-like coherent nonlinear structures [17]. 
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Figure Captions 
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Figure 1. Parameters T and a for the integrability of the coupled mode equations (14) 
and (15) with negative group dispersion. (3 = corresponds to the case of cold plasmas whereas 
large j3 corresponds to the weakly magnetized plasmas. For v = 6, parameters M and A can 
independently take arbitrary values for integrability. 



Figure 2. Parameters A and a for the integrability of the coupled mode equations (14) 
and (15) with negative group dispersion. (5 — corresponds to the case of cold plasmas whereas 
large j3 corresponds to the weakly magnetized plasmas. For v — 1 or 16, the corresponding 
values of 7 = ^fvY are to be obtained from Figure 1. Note that unlike the v = 6 case, both the 
parameters M and A are not arbitrary for v = 1 or 16, but are related by Eq. (22). 



